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INTRODUCTION 


In  this  report  we  summarize  the  results  of  exiK^riments  in  which  I.g  waveforms  arc  >y:j' !  ■  ' 
citlier  by  dynamic  raytracing  and  or  superptrsilion  of  Gaussian  beams  (Cerveny.  IbK'.  i  >  - 
dimensionally  varying  crustal  waveguides.  In  both  approaches,  l,g  is  rr-presenied  as  a  -  .c 
inulti()Iy  critically  rellecteri  SrnS  wav(-s  (Figure  I ).  Topography  of  i  lie  Molio  is  < orrelaled  wn  n  c  . .  i 
or  known  surface  topography  under  the  assumption  of  Airy  isosiasy.  The  starting  as'iinijc  , 
arc  essentially  the  same  as  those  use<l  by  K<'nneli  (1986)  in  id(‘ntifying  gross  eir.rts  of  m 
thickness  on  Lg  by  plotting  the  bounce  points  of  critically  reliectf'd  S  waves  in  the  crii't.  11":.  . 
iiowever,  waveforms  are  calculated  and  the  effects  of  geometric  sprearling  are  inchukvl.  Deta  ' . 
modeling  approach,  including  source  functions  assumtxl,  are  described  tti  {(’orinier  and  iOdir,!  . ,  <  i 
1991). 


RESULTS 

Propagation  Perpendicular  to  a  Mountain  Range 

bg  is  modeled  by  summing  dynamically  traced  rays  in  a  laterally  varying  crustal  waveguid.e  sium- 
lating  a  hypothetical  mountain  range.  Crustal  thickness  contours  are  shown  in  map  view  iti  I'ienie 
la.  Certain  characteristics,  such  as  the  square  corners  of  the  contour  lines,  are  concession>  m  idf 
for  modeling  simplicity  and  arc  not  expected  to  greatly  influence  the  results.  Background  < nstai 
lliickness  is  taken  as  35  km,  which  represents  continental  cru.st  slightly  elevated  above  .-ea  It  \(  1. 
'I’he  thickness  of  the  mountain  root  is  .50  krn,  which  is  consistent  with  gravity  snrvey.s  in  xatitii;-. 
high  mountainous  regions  (c.g.  Press  and  Sitwer,  1982,  p.  437).  .\  physical  analogy  to  i!i!.>  moiie! 
may  be  found  in  the  South  American  Andes,  which  is  similarly  a  narrow  range  with  abruj)!  bont' 
on  both  sides. 

Tested  propagation  paths  are  shown  in  Figure  la  and  ray  synthetic  seismograms  in  Figr.re  l!i 
for  a  vertical  point  force.  The  first  trace  show.s  Lg  at  the  .same  offset  as  receiver  Hj  throi;*:;!:  ;i>i 
undisturbed  waveguide.  'lYaces  2  through  4  are  recordings  from  rect'iver  locations  Rj  through  R 
Visual  comparison  of  the  traces  gives  a  measure  of  attenuation  and  the  extensive  nurnerii  a!  d.^ta 
from  the  program  runs  explain  the  primary  cau.se  (mantle  transmi.ssion  and/or  defocusing).  l  lic 
seismogram  from  location  Rj  shows  that  early  Lg  arrivals  have  high  amplitude  but  later  arrivals  arc 
considerably  weaker.  The  effects  of  crustal  thinning  are  not  evident  at  this  location  and  the  weaker 
signal  is  explained  entirely  by  propagation  across  the  transition  to  thicker  crust.  The  ttK  tiiod  of 
attenuation  is  defocusing.  The  slower  ray  paths  reaching  location  Ri  are  reflected  in  aretts  of  .Molto 
curvature  and  as  a  result,  average  geometric  spreading  is  increa.sed  by  a  factor  of  2.  Since  Moho 
curvature  is  found  to  cause  attenuation  in  this  and  other  program  runs,  it  is  useful  to  consider 
an  alternative  transition  model;  a  planar  dipping  Moho.  The  planar  surface  will  not  enuM'  itie 
increase  in  geometric  spreading  noted  above,  but  will  result  in  attenuation  through  other  int  ans. 
Spficifically,  much  of  the  defocused  energy  that  reaches  the  receiver  through  curved  tran.sii  ioii-.  i- 
reflected  by  planar  transitions  in  directions  that  do  not  reach  the  receiver.  This  is  etj-sily  \  i.'-.Mii/ed, 
by  considering  rays  tubes  and  the  area  between  rays,  which  is  the  wavefront.  Dcfoctising  I.  ,>.i- 
larger  wavefronts  which  have  a  better  chance  of  reaching  the  receiver  than  do  the  smalku  \va\  ( fi  out  s 


reflected  by  planar  surfaces.  The  present  study  has  compared  the  methorls  and  found  coinpaiahl. 
levels  of  attenuation,  although  the  seismograms  are  somewhat  different  in  apix-aram-e, 

A  very  slight  increase  in  attenuation  is  evident  at  receiver  location  H2,  primarily  Ixsau.'^e  .lii 
rays  reaching  this  point  are  defocused  to  some  e.xient.  Attenuation  at  location  R3  is  (juiti  :t  hn 
stronger  than  at  Ri  and  R2.  ffere,  the  propagation  path  crosses  the  trarihition  to  thini..  i  ciu.m. 
which  diverts  much  of  the  energy  that  would  otherwise  reach  the  receiver.  Diversion  (xems  in 
three  ways.  First,  rays  travelling  at  post-critical  angles  relative  to  the  horizontal  Moho  hf'couK 
pre-critical  at  the  upward  dipping  transition,  and  energy  is  lost  to  mantle  transmission.  SccoiuHn. 
a  large  portion  of  the  remaining  reflected  energy  travels  at  pre-critical  angles,  leading  to  mamli 
transmission  beyond  the  transition.  The  third  method  of  diversion  is  hack-reflectwl  rays. 

Propagation  Parallel  to  a  Mountain  Range 

Figure  2  illustrates  the  results  of  tests  of  Lg  propagation  parallel  to  the  strike  of  a  inouiiin:;, 
range  for  a  vertical  point  force.  Sources  an<l  receivers  are  assumed  to  lie  within  a  monnlain  taitgi 
having  a  deep  crustal  root.  Seismogram  trace  1  shown  again  for  reference  |)ropagation  in  a  < m--! 
having  planar  free  surface  and  Moho.  Traces  2-4  correspond  to  receiver  locations  Rj,  Ra  anri  H  ,. 
Transition  zones  trending  parallel  to  the  propagation  path  form  a  secondary  waveguide  and  1  hi 
modeling  predicts  that  rays  will  be  trapped  within.  These  raypaths  are  weak,  however,  as  Moho 
curvature  causes  significant  attenuation  by  geometric  spreading.  The  seismogram  at  location  Uj. 
which  is  not  subject  to  transitions  perpendicular  to  propagation,  shows  energy  roughly  ecpial  to 
undisturbed  Ig.  Although  a  large  number  of  raypaths  with  lateral  reflections  we,re  traced,  the 
increased  geometric  spreading  leads  to  weak  arrivals.  The  traces  at  locations  R2  and  R.3  are  .s.x'n 
to  be  strongly  attenuated  by  the  crustal  thinning  transition. 

Rays  trapped  by  the  lateral  waveguide  undergo  conversion  betwetm  SV  and  SlI,  but  conversion^ 
of  this  type  are  not  noticeable  in  the  synthetic  seismograms.  It  is  possible?  that  because  the.se  rays 
tend  to  be  weak,  the  conversions  are  not  noticeable.  Another  possibility  is  tliat  Sll  to  SV  convi'rsion 
is  generally  canceled  by  SV  to  Sll  conversion  under  the  particular  conditions  of  Ihi.s  model.  In  1  In 
real  crust,  the  primary  method  of  SV  to  SH  conversion  is  scattering  by  small  scale  heterogeaieit  ii  >. 
which  can  not  be  considered  under  the  constraints  of  ray  theory. 

Propagation  Oblique  to  a  Mountain  Range 

Figure  3  illustrates  the  attenuation  of  Lg  for  paths  oblique  to  the  strike  of  a  mountain  range  lor 
a  vertical  point  force.  Sources  and  receivers  are  located  such  that  proi)agation  paths  an-  at  iX), 
70,  45  and  20  degrees  relative  to  the  waveguide  perturbations.  Receiver  distance  is  2000  kin  in 
all  ca,ses.  Strong  attenuation  is  evident  over  each  propagation  path,  but  the  individual  wavefornis 
vary  considerably.  For  paths  terminating  at  receivers  Ri  and  R2,  rays  encounter  slot'p  Ihicknnss 
transitions  and  many  rays  are  entirely  transmitted  to  the  mantle.  Along  the  more  ohlicjnc  p.itli,; 
leading  to  R3  and  R4,  the  Moho  dips  less  steeply  and  as  a  result,  the  munlxir  of  two-point  ray  pal  hs 
increases.  It  appears  that  Lg  propagates  slightly  more  efficiently  at  45  degrees  than  at  the  oilx  t 
angles  tested,  but  the  diffonince  is  insignificant  compared  to  the  level  of  attenuation. 
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Width  of  Crustal  Transition 

In  crustal  models  shown  in  figure  la,  4a  and  41),  the  width  of  the  transition  region  thick 

and  thin  crust  is  variotl.  The  sourco-n'ceiver  paths  are  2000  krn  in  length  and  pf-rpendit  n:!!  n 
the  strike  of  the  transition  regions.  Syntlx'tic  traces  2  and  3  in  l-'ignre  Ic,  which  conipiU(  a  .'lio 
krn  width  of  thickened  crust  with  a  fiOO  ktn  width,  lx)th  show  weak  I,g.  The  similarity  lx  i  v-rt  ti 
traces  suggests  that  the  width  of  a  mountainous  region  has  very  li.,Ie  effect  on  overall  attenttation. 
'lYace  4  shows  Lg  after  propagation  through  transition  regions  200  krtt  wide  (compare<i  to  Mid  km 
wide  in  trace  2)  and  the  difference  is  readily  apparent.  The  widemsi  transition  width  re.Nuh.'.  in  Lg 
amplitude  analogous  to  that  in  the  undisturbed  waveguide,  which  im[)lie.s  that  the  wirliii  of  ( r'.i>t;i! 
transition  regions  and  the  corresponding  change  in  overall  Moho  dip  angle  impose  the  primarx 
control  over  Lg  attenuation. 

Propagation  Across  the  Tibet  Plateau 

The  map  shown  in  Figure  5a  gives  source  and  receiver  locations  tised  in  the  I,g  analysis  of  ILi/.iikni 
ot  al.  (1977).  Shading  corresponds  to  regions  of  high  mountains.  Central  Asia  wa.«-  chose!:  to 
test  the  methods  of  the  present  study  becaus«>  Lg  is  known  to  be  highly  atlonnaUsi  here  and  good 
specific  examples  arc  given  in  published  .seismograms.  A  study  of  ftKal  mechanisms  in  the  'ril)etan 
Plateau  (Molnar  and  Chen,  1983)  found  these  earthquakes  to  occur  at  depths  of  5  to  10  km  -.vith 
combined  strike-slip  and  normal  faulting  and  T  axis  primarily  <'ast-west.  No  earthciuakc.s  wer“ 
located  below  10  to  15  km,  suggesting  the  crtist  i.s  essentially  a.sei.smic  Ix'low  the.se  deotl;-.,  I  hc 
fault  solutions  indicate  the  Tibetan  Plateau  is  undergoing  ca.st-w(>st  extension,  much  like  the  Ihiviii 
and  Range  province  of  the  Western  United  States. 

Crustal  thickness  has  been  hand  digitiztxl  and  the  results,  in  the  form  of  a  Moho  dey»th  plot, 
are  shown  in  Figure  5b.  Samples  are  taken  at  a  .5  degree  (50  km)  sampling  of  latitude  and 
longitude.  (Topography  and  Moho  depth  at  10  kra  sampling  are  now  available  in  uncla.ssjfiod 
databases  established  by  Fielding  et  al.,  1992.)  The  method  of  determining  crustal  thickness  i.'  a 
combination  of  published  survey  results  and  inference.  Choudary,  1975  estimates  crustal  thicknes.s 
of  70-72  km  in  the  Central  Himalayas  from  Bouger  gravity  anomalies.  Bird  and  Tok.soz,  1975  stud> 
the  velocity  of  20-80  second  Rayleigh  waves  and  conclude  that  the  crust  is  75  km  thick  within 
the  Til>etan  Plateau.  Published  data  on  crustal  thickness  within  the  Tarim  Basin  and  Tien  Shiin 
Mountains  were  not  found,  so  thickness  has  been  estimated  from  average  elevation.  Flexation 
varies  from  500  to  1000  meters  within  the  Tarim  Basin  and  crustal  thickness  is  c.stimate(]  m  35 
km.  Average  Tien  Shan  elevation  is  on  the  order  of  3000  meters,  leading  to  a  thickness  estimate 
of  50  km.  Estimation  of  crustal  thickness  combined  with  the  large  digital  sampling  intei  vals  lead 
to  a  very  rough  mtKlol  of  the  crust,  but  the  tran.sitions,  which  seem  to  most  strongly  effec  i  Lg 
propagation,  arc  repr<'.s('ii(.(xl  as  clostjly  its  the  iruKiel  will  allow. 

Figure  5a  shows  propagation  paths  froni  tlu!  two  earthquakes  chosen  for  modeling  by  the  ray 
trifithod.  Double  couple  radiation  patterns  have  Ix'cn  includt'd.  Both  paths  are  subject  to  cni-,!;!! 
thinning  at  the  Tarim  Basin  then  crustal  thickening  at  the  Tien  ,Shan  mountains  before  the  waves 
arc  rrreorded  at  Soviet  station  Talgar  (TLC).  Since  the  paths  art'  similar,  transition  ste('()ne.vs  will 
1)0  approximately  the  same.  The  main  difference  between  paths  is  that  the  the  more  westerly  path 
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from  the  earthquake  near  87"K  longitude  travels  about  half  as  far  within  the  'I’ihwtan  Idateau  a> 
the  path  from  the  earthquake  near  93“  longitude.  Kigures  6a, b  compare  the  observcHi  and  s\  ni  liei  ii 
seismogram  for  the  earthquake  (mb--5.1)  having  the  more  westerly  path.  Lg  appears  to  be  sironeK 
transmitted,  having  a  high  amplitude  relative  to  Pn.  The  synthetic  agr(,>es  well  with  details  of  th. 
observed  Lg  coda.  Figures  7a, b  compare  seismograms  recorder!  from  earthquake  (mb  6.1)  having 
the  more  the  easterly  path.  Lg  appears  to  be  very  weak  in  these  recordings,  having  a  low  arnplinuie 
relative  to  Pn.  The  synthetic  seismogram  has  high  amplitude  pulses  not  .serm  in  the  observed 
seismogram. 

Certain  additions  to  the  crustal  model  could  significantly  improv'e  the  match  between  t  he  .syn¬ 
thetics  and  the  recordings,  hirst,  inclusion  of  the  low  Q  believed  to  exist  in  the  deep  crust  of 
the  Tibetan  Plateau  (Bird  and  Tok.soz,  1975)  would  improve  the  agreement  between  obsr'rved  an<i 
predicted  Lg  transmi.ssion  Each  ray  incident  at  th(>  Moho  In-ncath  tlu*  Tibetan  Plateau  would 
be  subject  to  high  intrinsic  attr-nuation,  establishing  a  direct  relationship  Ixttween  travel  dislatu  i 
within  the  plateau  and  attenuation  sugge.sted  by  the  observations.  Scattering  attenuation  by  tlie 
3-D  heterogeneities  associated  with  the  uplifted  region  would  tend  to  reinforce  this  relationsiiip 
Details  of  the  shallow  structure  of  the  Tarim  Basin  have  been  omittrxl  from  the  crustal  model. 
The  weak  path  ,  for  which  poor  agreement  exists  between  observed  and  pre(iicted  Lg  transmission, 
crosses  the  Tarim  Basin.  Baumgardt  (1991)  found  that  sedimentary  basin  structure  is  at  !<  a>t 
as  important,  if  not  more  important,  than  Moho  topography  in  controlling  Lg  propagation.  Low 
.seismic  velocity  within  the  sediments  can  form  a  waveguide  trapping  a  significant  portion  of  tin 
rays  that  would  otherwise  contribute  to  Lg,  the  slopes  of  the  basin  on  exiting  side  rays  acting  to 
defocus  subsequent  SmS  rays  (Figure  8).  Rays  reverberating  within  the  sedimentary  basin  would 
also  suffer  high  levels  of  scattering  attenuation  within  the  sediments.  A  third  po.ssibiiity  for  the  dis- 
crepency  between  the  observed  and  predicted  Lg  waveforms  in  the  weak  pat  h  exampk'  is  suggesii'd 
by  comparing  the  topographically  inferred  Moho  with  that  determined  by  the  Soviet  Deep  Seismu 
Sounding  Project  (DSS).  In  the  region  of  the  Tibet  l'*lateau,  the  transition  zone  betwwn  uonnal 
and  thickened  crust  is  quite  different  in  form  in  DSS  sections  than  that  which  would  be  predielixl 
by  simple  isostasy  (Fielding,  personal  communication).  One  or  more  of  the.se  modifications  of  ilv 
crustal  model  may  lead  to  synthetics  that  more  closely  match  the  attimiiation  levels  .seen  in  lb. 
recorded  seismograms. 

Scattering  Within  the  Crustal  Waveguide 

The  modeling  described  in  the  previous  suhsections  does  a  remarkably  ginid  job  in  pn'dii  iing 
the  blocking  of  Lg  by  major  crustal  transitions  and  is  often  .-ucce.ssful  in  {iccouriting  for  the  nio^i 
significant  details  of  the  Lg  coda.  It  is  clear,  however,  that  observed  Lg  codas  will  lx?  always  coiilain 
greater  complexity  than  synthetic  codas  generated  in  models  having  a  honiogiuieous  crust.  ^\V■  haw 
corulucted  several  experiments  to  include  the  effects  of  fine  .scale  three-dirni'iisional  lu'l.erogeneii  . 
within  a  crustal  waveguide  of  variable  thickne.ss.  To  achieve  this  goal,  wc  added  dynamically  t  rai  <  .i 
rays  to  point  scatterers  within  a  crustal  waveguide  of  variable  thickness  (I-’iguro  9).  A  staiKiard 
Born  approximation  (e.g.,  Wu  anil  Aki,  1985)  is  usisl  to  calculate  the  radiation  pattern  of  the  point 
scatterers.  The  results  shown  in  Figure  9  were  calculatwl  for  a  single  point  heterogonoi*v  having 
positive  5  per  cent  contrasts  in  P  velocity,  S  velocity,  and  density  from  the  backgronnd  nuxliuni. 
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1(.  is  dear  that  the  existence  of  even  a  single  jxwnt  heterogeneity  (jiiickly  acts  to  iill  ui>  i!u'  <  ...ia 
between  multiply  reflected  S  wavtiS. 

The  standard  Born  approximation  use<l  in  this  demonstration  does  not  conserve  energy  arid  d  ><  s 
a  poor  job  in  recreating  scattered  codas  predicKxl  by  finite  diflerenoe  solutions  of  the  elasti<'  cosuu  ion 
of  motion  in  heterogeneous  mt*dia  (Coates  and  Charrette,  1992).  Improved  agreement  with  finio 
difference  simulations  can  be  obtainetl  by  using  the  generalized  Born  approximation  (Coao  <  aiKi 
Chapman,  1991).  The  generalized  Born  ap[)roximation,  in  which  .scattering  occurs  In-  a  psof  cs' 
of  frequency  dependent  reflection  by  regions  of  high  gradient,  is  applicable  to  models  in  w  hif  h 
velocities  and  density  are  parameterized  by  continuous  functions  of  space.  Such  a  parameK  ri/c.'inri 
also  greatly  simplifies  the  description  of  the  rays  contributing  regional  .seismograuis. 

Importance  of  Scattering  in  Different  Depth  Regions  of  the  Crust 

Lg  rays  have  been  dynamically  traced  in  moriels  of  the  crust  in  wdiich  the  Moho  has  been  rep:  e>e’ 
by  a  high  gradient  transition  (Figure  10).  The  SmS  rays  contributing  to  the  l,g  cod, a  are  uuuinu 
rays  in  the  high  gradient  region.  Geometric  spreading  factors  have  Ixen  tracked  along  ea(  h  rav 
demonstrating  that  the  spreading  factors  of  individual  rays  vanish  at  the  turning  points  v.nhii. 
the  Moho  transition.  It  is  well  known  that  such  a  high  gradient  zone  will  form  a  caustic  stirface 
within  the  high  gradient  zene,  to  which  multiple  turning  rays  will  lx-  tangent.  A  reversed,  tra'.e! 
time  branch  and  triplication  will  also  be  observed  for  rays  traced  at  a  sequence  of  increasing  i.ike- 
off  angles.  Scattering  of  Lg  rays  near  the  caustic  surface,  where  amplitudes  will  be  high  due  to 
the  vanishing  of  the  geometric  spreading  factor,  will  be  greatly  enhanced.  This  situation  is  very 
analogous  to  the  enhanced  amplitude  of  scattered  PKIKP  precursors  due  to  heterogeneity  at  the 
core  mantle  boundary  (e.g.,  Bataille  and  Finite,  1988).  For  PKIKP  precursors,  the  presence  of 
a  caustic  surface  in  the  outer  core  is  siifncient  to  make  even  a  very  slight  heterogeneit.\’  rx  .u  the 
core-mantle  boundary  account  for  the  ob.srTved  amplitudes  of  the  scattered  waves,  l.ikewisr.  a 
caustic  surface  within  a  Moho  transition  zone  will  capable  of  amplifying  the  effects  of  veiy  .-'ight 
!i(;terogeneity  within  the  Moho  transition. 

In  predicting  regions  in  which  enhanced  scattering  may  oerur  it  is  also  useful  to  consiriet  tin 
intersections  of  SmS  rays  traced  from  the  receiver  with  SmS  rays  that  originate  at  the  sourtc.  In 
two-dimensional  scattering,  the  existence  of  a  heterogeneity  at  such  crossing  points  will  i  iakc  it 
possible  for  SraS  to  SmS  scattering  to  occur,  filling  up  the  coda  between  the  individual  SmS  arrivals 
making  up  the  Lg  coda.  The  greatest  density  of  such  crossing  points  occurs  near  the  Moho  t  ransitiou 
and  near  the  free  surface.  An  important  .sourct;  of  the  late  Lg  coda  will  be  three-dimensional  SmS 
to  StnS  scattering  occurring  out  of  the  vertical  plane  containing  iIk'  source  and  receiver. 

CONCLUSIONS  AND  RECOMMENDATIONS 
Modeling  Results 

(1)  Moho  Topography  and  Sedimentary  Basin.s  are  the  two  mosi,  imttortant  elements  stnx  t  hi  - 
responsible  for  Lg  blockage. 


(2)  Slof)es  of  major  crustal  discontiruiitics  and  transitions  (Molio  and  basin-sc<iiiru‘nt  inlet  fuel  ei. 
the  principal  features  controlling  the  propagation  elficieiu  y  of  l>g. 

(3)  There  is  only  a  weak  depr-ndence  on  the  azimuth  of  the  propagation  path  with  respe<  t  te  iIh 
strike  of  blocking  features. 

(4)  Distributed  heterogeneities  arul  small  srale  to[K)graphy  of  crustal  discont intiities  in  'he  <  iit^! 
can  act  to  fill  up  and  lengthen  tlie  Lg  r-orla.  Ih'lerogeneities  within  the  mist  primarily  aione  il,. 
source-receiver  azimuth  fill  u[)  the  early  coda  betwet'ii  individual  S  waves  lra[)[)ed  in  the  <  ru'.iai 
waveguide.  Off  azimuth  heterogeneities  within  the  crust  are  a  significant  component  of  tlie  late  l  .j.; 
coda  following  the  last  arriving  .\loho  reflected  S  wave. 


Lessons  from  Ray  Based  Approaches  to  Regional  Seismogram  Synthesis 

(1)  Utility  of  a  Continuous  .Model  Parameterization:  The  complexity  of  ray  code  descriptions  in 
models  described  by  first  ordtT  discontinuities  argues  for  motlels  paramet('rize<i  by  functions  din 
tinuous  in  space. 

(2)  'freatment  of  Reflections  and  Conversions  in  Regions  of  Strong  Cradient;  A  geiKTaiizeil  Horn 
approximation  can  be  used  to  calculate,  the  frequency  dependent  scattering  by  regions  having  st  ronu 
gradient,  including  the  Moho  transition  zone  resulting  from  the  parameterization  describt'd  above. 

(3)  The  Most  Important  Regions  of  Scattering  for  the  Lg  Phase  and  Why; 

(a)  Moho  transition  b<t:ause  a  canstic  snrfa<-e  is  ineludfsi  in  this  region  for  many  of  tlie  S  w.im  ,-. 
('omi)osing  the  Lg  phase.  In  addition  there  is  a  high  d<>nsity  of  S  ray  inter. sections  in  this  region 
permitting  scattering  from  one  type  of  multiply  reflected  S  wave  composing  the  Lg  phase'  to  anoi  her 
type. 

(b)  The  free  surface  region  because  there  is  i  high  density  of  S  ray  intersections  in  thi.s  ic-gion 
permitting  scattering  from  a  source  ray  consisting  of  multiple  Moho  nrflocbHl  S  waves  to  a  reeeivet 
ray  consisting  of  multiply  reflecUxl  S  waves.  In  addition,  this  region  should  bi'  important  ba.seil  ou 
array  data  suggesting  a  higher  magnitude  of  velocity  and  density  nuctuations  in  the  nppt'r  2  km  ol 
the  crust. 


(i 


0 


500  km 


Figure  i  (a)  Above:  ()ropaganon  paths  norma!  to  strike  of  mountain  range.  This  map  view  shows 
the  crustal  thickness  contours  of  a  narrow  hypoihelicai  mountain  bell,  along  with  the  Lg  propaga¬ 
tion  path  and  receiver  locations.  Ri  through  Ra  offsets  are  2200,  2300,  and  2400  km  respectively,  an 
unperturbed  waveguide,  (b)  Below:  ray  synthetic  seismograms  for  paths  to  receivers  Ri  through 
R3.  'lYace  1  is  Lg  modeled  at  the  Ri  offset  in  an  unperturbed  waveguide.  Transverse  ground 
motion  is  not  shown  because  SV  to  SH  conversion  does  not  occur  when  the  propagation  path  is 
perpendicular  to  wavt'guido  variations. 


1  Minute 


Trace  1  =  Non -attenuated  Lg 
Trace  2  =  Receiver  R, 

Trace  3  =  Receiver  R, 

Trace  4  =  Receiver  R, 
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0  500  km 


Figure  2 

(a)  Above:  propagation  paths  parallel  to  the  strike  of  the  mountain  range.  Offsets 
are  2200,  2300,  and  2400  km.  to  receivers  within  the  mountain  range,  (b)  Below:  ray  sjmthetic 
seismograms.  An  increased  number  of  ray  paths  reach  receiver  Rj  relative  to  receivers  R2  and  R3 
resulting  from  lateral  reflections  within  the  mountain  belt.  Coda  duration  is  sharply  limited  at  the 
Rj  and  R3.  V'ery  weak  motion  can  be  seen  on  the  transverse  traces. 


1  Minute 


Trace  1  =  Non-atienuated  Lg 
Trace  2  =  Receiver  R, 

Trace  3  =  Receiver  R,  8 

Trace  4  =  Receiver  R, 


3  5  km 
50  km 


I  I  -  I _ l _ \ _ j 

0  500  km 


Figure  4  q{  width  of  crustal  thickening  and  thinning.  The  crust  increases  and  decreases  in 

thickness  over  transition  regions  100  km  (above)  and  200  km  (middle),  ilay  synthetic  Lg  wavetrains 
are  shown  below  for  an  offset  of  2000  km  and  propagation  perpendicular  to  the  strike  of  the  mountain 
range.  Weak  Lg  in  traces  2  and  3  compared  to  strong  Lg  in  trace  4  indicates  that  the  width  of  the 
transition  region  affects  Lg  efficiency  more  than  the  overall  width  of  the  mountain  root. 
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1  Minute 


Trace  1  =  Non-attenuated  La 
Trace  2  =  Receiver  R,  Crust  Model  la 
Trace  3  -  Receiver  R^  Crust  Model  4a 
Trace  4  =  Receiver  R,  Crust  Model  4b 
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Degrees  East 


“  Above:  (a)  map  of  Central  Asia  showing  the  propagation  paths  of  two  earthquakes 

modeled  by  ray  methods.  Shading  highlights  high  mountainous  regions  of  thickened  crust.  The 
left  path  gives  a  detectable  Lg  signal  and  is  referred  to  as  the  strong  path.  Lg  from  the  right 
earthquake  is  very  weak.  Below:  (b)  Central  Asian  Moho  depth  plot  obtained  by  digitization  of 
crustal  thickness  values.  Not  to  scale. 


Figure  6  Above:  (a)  Conirai  Asia  strong  path  .stnsmograms.  The  seismograms  recorded  from  the 
loft  earthquake  of  figure  5a  show  clear  Lg  energy  m  both  the  top  trace  i  instrument  response  peaked 
between  0.5  to  1  Hz.)  and  bottom  trace  (instrument  response  peaked  between  0.2  to  0.5  Hz.). 
Middle:  (b)  scaled  strong  path  synthetic.  The  synthetic  seismogram  is  scaled  to  the  dimensions  of 
the  actual  seismogram  and  shows  good  recreation  of  Lg.  Below  (c):  detailed  strong  path  synthetic. 
These  synthetic  seismograms  are  not  scaled  to  the  recordings  and  allow  determination  of  the  level 
of  attenuation.  Ray  modeling  through  a  homogeneous  crust  predicts  significant  attenuation  over 
this  path. 


10  Sec 


Trace  1  =  Non-attenuated  Lg 
Trace  2  =  Strong  path  synthetic 
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Figure  7  Above:  (a)  Central  Asia  weak  path  seismograms.  The  seismograms  recorded  from  the 
right  earthquake  of  figure  5a  show  very  weak  Lg,  as  can  be  seen  by  comparing  Lg  amplitude  to 
Pn.  Middle:  (b)  scaled  weak  path  synthetic.  The  synthetic  seismogram  of  the  weak  path  predicts 
energetic  pulses  that  are  not  present  in  the  actual  recordings.  Below:  (c)  detailed  weak  path 
synthetic.  Lg  is  clearly  attenuated  over  this  path  but,  as  seen  in  figure  9b,  the  predicted  level  of 
attenuation  is  insufficient  to  explain  the  observations. 


ta  Sac 


Trace  1  =  Non*attenuated  Lg 
Trace  2  =  Weak  path  synthetic 
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3da i n  Mode  i 


000. 


Lg  ray  trace 


'I'Bniiig  a  deep  sedimcniary 

basin  (Z  and  Xco-ordinaies  are  in  kilometers.)  Bottom:  ray 

rtM'ih/  waveguide  (Lg).  .Mole 

^ihe  effect  of  the  basin  will  be  to  attenuate  Lg  at  stations 
1400  km  and  greater  range. 
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Figure  9.  Effect  of  scattering  within  a  crustal  waveguide  of  variable  thickness.  Left 
column:  no  scattering  within  the  crustal  waveguide.  Right  column:  scattering  by  a  single 
scattercr  within  the  crustal  waveguide.  (In  the  synthetics  shown  in  the  right  column,  all  of 
the  ray  paths  shown  at  top  left  are  included  plus  the  scattering  paths  shown  at  top  right.) 
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Figure  10.  Intersections  of  rays  traced  from  the  receiver  with  rays  traced  from  the  .source 
are  densest  near  the  surface  and  near  the  SmS  turning  points.  Scaiterers  at  these  locaiion.s 
will  contribute  to  the  complexity  of  Lg  coda.  Scattering  near  the  Moho  transition  will  be 
particularly  important  becau.se  the  presensce  of  a  caustic  in  this  region  will  greatly  amplify 
the  effects  of  small  perturbations  in  elastic  moduli  and  density. 
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